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(57) Abstract: The invention concerns an adapted filter (5) producing observations of a reference vector (r) of W successive samples 
of a signal formed by adapted filtering of a signal (Y) received through a radio propagation channel with a pilot spread code (Cr). 
A processing is performed in said adapted filter output to estimate the parameters representing the propagation channel, including 
specific vectors (Vy) of a self-correlating matrix of the reference vector and the specific associated values (X^). Said parameters can 
be used to perform coherent demodulation of the signals received through other logic channels sharing the same physical propagation 
channel. 



(57) Abregd: Un filtre adaptd (5) produit des observations d’un vecteur de reference (r) de W dchantillons successifs d’un signal 
forme par filtrage adapts d’un signal (Y) regu selon un canal de propagation radio avec un code d’dtalement pilote (Cr). Un traitement 
est effectud en sortie de ce filtre adaptd pour estimer des paramfctres reprdsentant le canal de propagation, incluant des vecteurs 
propres (v) d’une matrice d’ autocorrelation du vecteur de reference et les valeurs propres associees (X k ). Ces param£tres peuvent 
£tre utilises pour effectuer une demodulation coherente de signaux regus suivant d* autres canaux logiques partagearit le meme canal 
physique de propagation. 
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DISPOSITIFS D’ESTIMATION DE CANAL RADIO ET DE DEMODULATION 

La presente invention concerne les techniques d’estimation des canaux 
radio, ainsi qu’un demodulateur utilisant une telle technique. 

L’invention trouve application dans le domaine des 
5 radiocommunications a acces multiple a repartition par codes (CDMA, « Code- 
Division Multiple Access »), notamment dans les cas ou le temps d'un symbole 
d’information est sensiblement plus grand que la memoire du canal (c'est-a- 
dire qu’on peut negliger (’interference inter-symboles). 

Dans un systeme CDMA, les symboles transmis, binaires (±1) ou 
io quaternaires (±1 ±j), sont multiplies par des codes d’etalement composes 
d’echantillons, appeles « chips », dont la cadence est superieure a celle des 
symboles. Des codes d’etalement orthogonaux ou quasi-orthogonaux sont 
alloues a differents canaux logiques partageant la meme frequence porteuse, 
afin de permettre a chaque recepteur de detecter la sequence de symboles qui 
15 lui est destinee, en multipliant le signal regu par le code d’etalement 
correspondent. 

Un avantage du systeme CDMA est qu’il permet d’exploiter des trajets 
multiples de propagation radio, en utilisant ce qu’on appelle un recepteur en 
rateau (« rake receiver »). 

20 Le recepteur en rateau effectue une demodulation coherente fondee 

sur une approximation de la reponse impulsionnelle du canal de propagation 
radio par une serie de pics, chaque pic apparaissant avec un retard 
correspondant au temps de propagation le long d’un trajet particulier et ayant 
une amplitude complexe correspondant a I’attenuation et au dephasage du 
25 signal le long de ce trajet (realisation instantanee du fading). En analysant 
plusieurs trajets de reception, c’est-a-dire en echantillonnant plusieurs fois la 
sortie d'un filtre adapte au code d’etalement du canal, avec des retards 
correspondant respectivement a ces trajets, le recepteur en rateau obtient des 
estimations multiples des symboles transmis, qui sont combinees pour obtenir 
30 un gain en diversity. La combinaison peut notamment etre effectuee selon la 
methode dite MRC (« Maximum Ratio Combining »), qui pondere les 
differentes estimations en fonction des amplitudes complexes observees pour 
les differents trajets. 

Afin de permettre cette demodulation coherente, un canal pilote peut 
35 etre prevu pour I’estimation de la reponse impulsionnelle sous forme d'une 
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succession de pics. La reponse impulsionneile est estimee au moyen d’un filtre 
adapte a un code d’etalement pilote avec lequel I’emetteur module une 
sequence de symboles connus, par exemple des symboies a 1. Les positions 
des maxima de la sortie de ce filtre adapte donnent les retards utilises dans les 
5 doigts du recepteur en rateau, et les amplitudes complexes associees 
correspondent aux valeurs de ces maxima. Lorsque plusieurs trajets distincts 
ont des retards tres proches, par exemple avec des ecarts inferieurs a la duree 
d’un chip, la recherche des maxima conduit generalement a ne prendre en 
compte qu’un seul de ces retards. 

io Un but de la presente invention est d’exploiter une representation du 

canal de propagation radio qui soit moins approximative que dans un recepteur 
en rateau traditionnel, notamment dans le cas ou il existe plusieurs trajets de 
propagation ayant des retards proches. 

Selon I’invention, il est propose un dispositif d'estimation de 
15 parametres representant un canal de propagation radio, comprenant un filtre 
adapte pour produire des observations d’un vecteur de reference de W 
echantillons successes d’un signal forme par filtrage adapte d’un signal regu 
selon le canal de propagation avec un code d’etalement pilote, et des moyens 
d'estimation de parametres representant le canal de propagation, incluant des 
20 vecteurs propres d’une matrice d’autocorrelation du vecteur de reference. Les 
parametres representant le canal de propagation peuvent en outre comprendre 
des valeurs propres de la matrice d’autocorrelation du vecteur de reference, 
respectivement associees auxdits vecteurs propres avec lesquels elles sont 
estimees. 

25 Le dispositif emploie une description du canal de propagation sous 

forme de ses elements propres, ce qui procure une estimation plus riche que le 
modele de la ligne a retard (retards et amplitudes) pris en compte dans le 
recepteur en rateau traditionnel. Ainsi, chaque trajet de propagation peut etre 
pris en compte dans la reponse du canal au moyen d’une forme d’onde 
30 correspondant a un vecteur propre de la matrice d’autocorrelation et non 
simplement par un pic ponctuel associe a un retard. 

Dans un mode de realisation prefere, I’estimation des vecteurs propres 
de la matrice d’autocorrelation comprend les etapes suivantes : 

- extraire du vecteur de reference M sous-vecteurs de reference composes 
35 chacun de S echantillons successifs dudit signal forme par filtrage 

adapte, M et S etant des nombres plus grands que 1, les M sous- 
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vecteurs de reference correspondent en moyenne a M plages 
energetiques du vecteur de reference ; 

estimer les matrices d autocorrelation respectives des M sous-vecteurs 
de reference ; 

® ” determiner des sous-vecteurs propres associes chacun a une valeur 

propre et a I’un des M sous-vecteurs de reference, chaque sous-vecteur 
propre etant, pour la valeur propre associee, un vecteur propre de la 
matrice d autocorrelation estimee du sous-vecteur de reference associe ; 

- former chaque vecteur propre estime de la matrice d’autocorrelation du 
10 vecteur de reference en plagant un des sous-vecteurs propres aux 

positions des S echantillons du sous-vecteur de reference associe et des 
echantillons nuls aux W— S autres positions. 

Ceci permet de realiser un compromis entre la finesse d'estimation du 
canal et la complexity des calculs requis, qui sont toutes deux des fonctions 
15 croissantes de S, ainsi que de M. Plus la capacity de calcul disponible est 
grande, plus les sous-vecteurs peuvent etre de taille importante, ce qui affine 
I’approximation de la reponse du canal. 

Le dispositif ci-dessus peut etre utilise chaque fois qu’on a besoin de 
caracteriser un canal de propagation, en particular pour effectuer une 
20 demodulation. 

Un autre aspect de la presente invention se rapporte ainsi a un 
demodulateur pour une reception selon un canal de propagation radio sur 
lequel sont definis plusieurs canaux logiques multiplexes incluant un canal 
pilote et au moins un canal de donnees, comprenant des moyens d'obtention 
25 de plusieurs observations d'un premier vecteur de W echantillons successifs 
d un premier signal forme par filtrage adapte d’un signal regu selon le canal de 
propagation au moyen d’un premier code d’etalement relatif au canal pilote, 
des moyens d’obtention d’un second vecteur de W echantillons successifs d’un 
second signal forme par filtrage adapte dudit signal regu au moyen d’un second 
30 code d’etalement relatif au canal de donnees, des moyens d’estimation de 
parametres representant le canal de propagation, incluant des vecteurs 
propres d une matrice d’autocorrelation du premier vecteur, et des moyens 
d'estimation de symboles pour estimer au moins un symbole d'information 
porte par le second signal sur la base des premier et second vecteurs et 
35 desdits parametres representant le canal de propagation. Ce demodulateur 
incorpore un dispositif d’estimation de canal tel qu’enonce precedemment. 
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Si le second signal porte des symboles d’information binaires modules 
selon une modulation binaire a deplacement de phase, les moyens 
d’estimation de symboles peuvent etre agences pour estimer au moins un 
symbole d’information binaire porte par le second signal par le signe de la 
5 partie reelle du nombre complexe : 




ou d est le nombre de vecteurs propres estimes de la matrice d’autocorrelation 
du premier vecteur, v k est un vecteur propre de la matrice d’autocorrelation du 

premier vecteur, associe a une valeur propre X k , N 0 est une estimation de la 
10 puissance du bruit sur le canal de propagation, r et x sont respectivement les 
premier et second vecteurs, (3 est un coefficient predetermine, (,) H designe le 
transpose conjugue, et ||.|| la norme d’un vecteur. 

Une vraisemblance dudit symbole d’information binaire estime peut 
etre determinee proportionnellement ci la valeur absolue de la partie reelle du 
15 nombre complexe y. 

Si le second signal porte des symboles d’information quaternaires 
modules selon une modulation quaternaire a deplacement de phase, les 
moyens d’estimation de symboles peuvent etre agences pour estimer au moins 
un symbole d’information quaternaire porte par le second signal par deux bits 
20 respectivement d^finis par les signes des parties reelle et imaginaire du 
nombre complexe y defini par la relation (1). Des vraisemblances respectives 
de ces deux bits peuvent en outre etre determinees proportionnellement aux 
valeurs absolues des parties reelles et imaginaires du nombre complexe y. 

Le demodulateur peut egalement etre adapte pour traiter un signal regu 
25 selon le canal de propagation sous forme de plusieurs signaux captes au 
moyen d'antennes differentes, afin de procurer une diversity spatiale. 

D'autres particularites et avantages de la presente invention 
apparaTtront dans la description ci-apres d’exemples de realisation non 
limitatifs, en reference aux dessins annexes, dans lesquels : 

30 - la figure 1 est un schema synoptique d’un exemple de recepteur de 

radiocommunication comportant un demodulateur selon I’invention ; 

- la figure 2 est un schema d’un module d’estimation de canal du 
demodulateur de la figure 1 ; 




WO 01/41382 



PCT/FROO/03306 



-5- 

- la figure 3 est un schema synoptique d’un module d’estimation de 
symboles du demodulateur de la figure 1 ; et 

- la figure 4 est un schema synoptique d'une unite de calcul du module de 
la figure 3. 

Le recepteur represente sur la figure 1 comporte un etage radio 1 qui 
re$oit le signal radio capte par I’antenne 2 et le convertit en bande de base. Le 
signal en bande de base est numerise par un convertisseur analogique- 
numerique 3, puis fourni a un filtre de reception 4. Le filtre 4 assure un filtrage 
adapte a la mise en forme des signaux par I’emetteur. II delivre un signal 
numerique Y a raison d’un echantillon complexe par chip des codes 
d’etalement. 

En presence de trajets multiples de propagation (au nombre de d), la 
reponse h(t) du canal allant jusqu’au filtre de reception 4 peut s’ecrire : 

d-l 

h(t)= £s k .a(t-x k ) (2) 

k=0 

Les t k represented les differents retards des echos, les s k sont les realisations 

instantanees du fading variant au cours du temps, et a(t) designe I’impulsion 
formant la modulation. On suppose que le critere de Shannon est respecte 
lorsque les signaux sont echantillonnes a la frequence chip 1/T, et que la 
reponse du canal est nulle pour t > (W-1)T. Dans ce cas, la reponse du canal 
peut etre representee par un vecteur a W composantes complexes : 

h(0) 

s= h(T > 

_h((W-1)T) 

Dans la relation (3), les sont les vecteurs propres de la matrice 

d’autocorrelation de h (soit B.v k = p k .v k avec R = E^).h H ), ou E(.) designe 
I’esperance mathematique), et les a k sont des realisations instantanees de 
variables aleatoires gaussiennes complexes de variances egales aux valeurs 

propres p k de la matrice d’autocorrelation (E^|a k | 2 j = |a k ). Cette decomposition 

de la matrice B donne des vecteurs propres orthogonaux (v k H v k . = 0 si k * k') 
et des valeurs propres ja k reelles positives ou nulles puisque la matrice est 
hermitienne. 



d-1 

k=0 



( 3 ) 
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Parmi les canaux logiques multiplexes sur le canal radio, on a un canal 
de controle sur lequel sont transmis des symboles pilotes et un ou plusieurs 
canaux de donnees sur lequel sont transmis des symboles d’information. On 
suppose que la puissance d’emission sur le canal de controle est (B fois 
5 superieure a la puissance d’emission d’un symbole d’information b sur un canal 
de donnees (cette derniere puissance sera consideree comme unitaire). Le 
coefficient 0 est connu a I’avance. Le signal est affecte d’un bruit blanc additif 
centre gaussien de variance N 0 . 

Un filtre adapte 5 filtre le signal regu Y conformement au code 
io d’etalement Cr du canal pilote. II en resulte pour chaque symbole pilote un 
vecteur r ayant W composantes complexes, donne par : 

d-1 

r = Vp£a k .v k +n (4) 

k=0 

avec E(n.n H )=N 0 .l , la matrice identite etant notee I . 

Un autre filtre adapte 6 filtre le signal regu Y conformement au code 
15 d’etalement Cx du canal de donnees. II en resulte pour chaque symbole 
inconnu b un vecteur x ayant W composantes complexes, donne par : 

d-i 

x = b.^a k .y k +n’ (5) 

k=0 

avec E(n'.n ,H )= N 0 .l et E(n.n' H )= 0. 



On note X = le vecteur gaussien de taille 2W constitue par la 

concatenation des signaux regus sur les deux canaux. Sa densite de 
probabilite p(X) peut s’ecrire : 



7^ w .det(K) 



P(X) = exp(-X H K -1 x) (6) 

7t 2W .det(K) v ~ ' 

ou K = E^C h x), et det(.) designe le determinant. 

Le but du demodulateur est de maximiser par rapport a b la fonction de 
vraisemblance donnee par la probabilite conditionnelle p(b|x). 

Comme les realisations du bruit et du fading sont independantes, on a : 
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K = 



PR b‘Vp.R 
b^/p.R R 



-7- 
d-1 



+ N 0 I = B. 



S( 2 ^k-^k-Wk) 



k=0 



B H + N 0 I 



(7) 



en utilisant la decomposition en elements propres de la matrice B, avec 



. Par application du lemme d’inversion 



2 


, 1 


Vk 


0 = b.l 


et w k = — p= . 

“ V2|v k 





B = 

matricielle, on obtient 



K -1 = 1 



Nr 



d-1 

!-Z 






( 8 ) 



k=o N 0 + 2^ k .|Bw k | 

En prenant le logarithme de la probability (6), et en notant que le 
determinant det(K) est commun a tous les symboles possibles et que 
2 

j|Bw k || =(l+p)/2, le demodulateur doit minimiser la fonction suivante par 
rapport a b : 



d-i 

4>(b) =£ 



2 Hk 



k e 0 N o+(i+PK 

a -' M K /||y k | 2 



X H Bw k 



k?o "o 



|vk-(#r + b.x)|‘ 



(9) 



Nn +(1+3)|ik 

Lorsque les symboles sont des bits modules selon une modulation 
binaire a deplacement de phase (BPSK, « binary phase shift keying ») ; 
I’estimation b d’un symbole b, qui minimise la fonction <|>(b), est alors donnee 
par le signe de la partie reelle du nombre complexe : 



d-i 



- kto'^o 



M-k / ||Yk 



•Kkx) 



(10) 



Nn +(1 + P)-M-k 

(b = sgn[Re(y)]), et la vraisemblance maximisee est proportionnelle a la valeur 
absolue de Re(y). Si un decodage a entrees souples (« soft input decoding ») 
est effectue en aval du demodulateur, la vraisemblance du bit estime b peut 
done etre mesuree par A = |Re(y)| . 

Lorsque les symboles sont quaternaires et modules selon une 
modulation quatemaire a deplacement de phase (QPSK, « quaternary phase 
shift keying »), I’estimation d’un symbole b = b 1 + j.b 2 , qui minimise la fonction 

4>(b), est formee par deux bits donnes par les signes respectifs des parties 
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reelle et imaginaire du nombre complexe y defini par la relation (10) 
(b = b-, + j.b 2 = sgn[Re(y)] + j.sgn[lm(y)]), et les vraisemblances de ces deux bits 
sont respectivement proportionnelles aux valeurs absolues des parties reelle et 
imaginaire de y. Si un decodage a entrees souples est effectue en aval du 
demodulateur, les vraisemblances des deux bits estimes b 1 ,b 2 peuvent done 
etre mesurees par A-) = |Re(y)| et A 2 = |lm(y)( . 



Lorsque le recepteur comporte deux antennes en entree de deux voies 
de traitement, chaque voie ayant deux filtres adaptes I'un au code Cr du canal 
pilote fournissant un vecteur r n , i^, et I'autre au code Cx du canal de donnees 

fournissant un vecteur x n , x^ la fonction a minimiser <J>(b) peut tenir compte de 



la diversity spatiale, son expression n’etant plus donnee par (9) mais par 

■*- /n n2 ^ _ /ii 1, 2 , ^ 



d -1 



«Mb)= X 



Hk/fekl 



k=ol N 0 +(1 + P)^ 



Yk- (VPC’+bx 1 ) 



4k/ Yk 



Nq + 0 + P) 4 k 



Yk H (VPC‘+bx') 



avec les notations : r' = 



fr n 

-n 






, x' = 



f x ^ 

-n 



K-dJ 



, R’ = 



K p *K 
p.K K 



( 11 ) 

, p etant le facteur de 



correlation entre les deux antennes (soit E(ak >n .a kd )= p.p k ). les vecteurs v k 



et v k etant les vecteurs propres de R‘ respectivement associes aux valeurs 
propres |i k + et p. k “ soit vj 

4k = ( 1 + |p|)4k et pk = (l -|p|Vk P° ur 0 < k < d 



1 


f Y k ] 


1 


f Yk ] 


V 2 ‘ 




, Vk = 1 

- k V 2 


l-PYkH 



Dans le cas de reppnses impulsionnelles de petite dimension W, il est 
aise d’extraire les vecteurs propres de R par une methode algebrique classique 
travaillant directement sur une estimation de la matrice, par exemple la 
methode dite « Power Iteration » (voir Golub & Loane : « Matrix Computation », 
editions Baltimore M.B.). 

Dans le cas de reponses impulsionnelles de plus grande dimension W, 
la complexite requise pour mettre a jour I’estimation de R rend cette methode 

irrealisable en pratique. On utilise alors une matrice R = E^.r H ] definie comme 
I’autocorrelation du vecteur r, liee a R par la relation : 
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R = p.R + N 0 .l 



( 12 ) 



On verifie immediatement que R et B ont les memes vecteurs propres 



v k . Les valeurs propres A. k de R sont de celles de R multiplies par (3 et 
decalees de N 0 , soit X k = p.p k + N 0 . Avec (10), on en deduit I’expression (1) du 

5 nombre complexe y servant a estimer les bits d’information. 

En general, la matrice R a une structure creuse et possede peu 
d’elements significatifs. Elle peut etre decomposee selon (’approximation 
suivante : 



00 ••• 0 
0 R° 

0 : 

'• 5 m_1 2 

o oo 



(13) 



10 



15 



c’est-a-dire sous forme de blocs de taille SxS, notes R' pour 0 < i < M, 
distribues le long de sa diagonale principale (1 < S « W). Les positions de ces 
M blocs le long de la diagonale de R correspondent aux zones les plus 
energetiques de la reponse impulsionnelle du canal. De meme, la matrice R 



peut etre approchee par : 

'N 0 .l 0 



R = 



0 R 



o 



Nn.l 



- M-1 

. R 0 

0 Nn.l 



(14) 



en notant R 1 = R 1 +N 0 .I I'observee bruitee de R 1 pour 0 < i < M. 



Chaque sous-matrice R 1 admet une decomposition en vecteurs 

propres, identique a celle de R‘, qui peut etre obtenue rapidement par une 
methode classique puisque sa dimension S est petite devant celle de R. En 
20 notant d' la dimension du sous-espace signal relativement a la sous-matrice 
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M-1 



R , avec ^d' = d, cette decomposition de R s’ecrit : 
i=0 ” 



d - 1 - r ■ 

^ j 



Si = I *-V 

j=0 



(15) 



De meme que precedemment, les valeurs propres sont decalees de la 
puissance du bruit. Ainsi les vecteurs propres v k de R sont : 

0 



v k = 



*1 



avec < 



fO < k < d- 1 
[0 < i < M, 0 < j < d 1 



(16) 



10 



15 



20 



et ses valeurs propres n k sont celles des M matrices R 1 , decalees de N 0 : 



Hk = 



O-k-Nq) ( x 'r N °) 

0 0 



avec < 



(0 < k < d-1 
0<i<M, 0<j<d' 



(17) 



Les sous-matrices R 1 sont localisees la ou la reponse impulsionnelle 

concentre ses maxima d’energie. Pour avoir la position des M sous-matrices, 
supposees de taille fixe SxS, il suffit done d’avoir une estimation du profil en 
puissance du canal et d’en extraire les M portions de taille S les plus 
energetiques. Ce profil peut etre simplement calcule a partir de la moyenne du 
module carre de chaque echantillon du vecteur r. Pour calculer les sous- 

matrices- d’autocorrelation R 1 , on determine les M sous-vecteurs H de 

dimension S correspondant aux signaux regus sur le canal pilote aux positions 
precedemment determinees, et on calcule apres un certain nombre 

• *H 

d’observations des moyennes des sous-matrices r'.r 1 . 

II se peut que certaines valeurs propres estimees ne soient pas 
significatives ou que Ton choisisse de ne retenir qu’un certain nombre fixe 
d'elements propres. On peut done, en fin de traitement, effectuer une selection 
des vecteurs propres associes aux plus puissantes valeurs propres. 

Le demodulateur de la figure 1 comprend un module d’estimation de 
canal 8, represente de fagon plus detaille sur la figure 2, qui calcule des 
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estimations des vecteurs propres v k et des valeurs propres X k de la matrice 
d'autocorrelation R du vecteur r delivre par le filtre adapte 5, ainsi qu’une 
estimation de la puissance du bruit N 0 . 

En reference a la figure 2, le module 8 comprend une unite 80 qui 
5 calcule des vecteurs de W composantes constitutes par les modules carres 
des composantes des observations successives du vecteur r delivre par le filtre 
adapte 5. Une unite 81 calcule une moyenne courante du vecteur delivre par 
I’unite 80, par exemple sur une fenetre rectangulaire ou exponentielle. (.’unite 
82 determine les M plages de S echantillons les plus energetiques du vecteur 
io moyenne, afin de definir les positions des sous-vecteurs r' de taille S (0 < i < M) 
au sein des vecteurs r. 

Une unite 83 extrait les M sous-vecteurs r* de chaque vecteur r, et en 
estime la matrice d’autocorrelation par un calcul de moyennes des produits 
terme a terme : 

-i ( i 

15 R = El r'.r J (19) 

Des unites de diagonalisation 84 extraient les valeurs propres xj des 

matrices R\ du moins celles de plus grands modules, et les vecteurs propres 
associes yi. Ce calcul peut etre effectue selon toute methode algebrique 

connue, par exemple « Power Iteration », qui ne pose pas de probleme de 
20 mise en oeuvre compte tenu de la reduction operee sur les dimensions des 
matrices a diagonaliser (SxS au lieu de WxW). L’unite 85 selectionne les d 
valeurs propres de plus grands modules X k parmi les valeurs propres X-j 
foumies par les unites 84, et les vecteurs propres correspondants v k qu’elle 
construit a partir des Vj* conformement a la relation (16). 

25 L’unite 85 peut egalement determiner la puissance du bruit N 0 , par 

exemple prise egale a la plus petite des valeurs propres X- D’autres methodes 
classiques d’estimation de N 0 peuvent etre utilisees. 

Revenant a la figure 1, les parametres v k , X k (0 < i < d) et N 0 estimes 

par le module 8 sont fournis au module 9 qui estime les symboles transmis a 
30 partir des vecteurs r et x delivres par les filtres adaptes 5, 6. 

La figure 3 illustre une structure possible du module 9 dans le cas 
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d’une modulation QPSK. Un additionneur 91 calcule la quantite y selon la 
relation (1) a partir des d contributions a sommer, calculees par des unites 
respectives 90. 

Une telle unite 90, relative a un trajet k (0 < k < d), est illustree par la 
5 figure 4. Les produits scalaires (r H .v k ) et (v k H .x) sont respectivement calcules 
en 95 et 96, puis multiplies entre eux par le multiplieur complexe 97. Le bloc 98 

effectue le calcul du coefficient ^ — — (on suppose ici que les unites 

(l+p)A k -(3.N 0 

84 ont delivre les vecteurs propres sous forme orthonormee, soit 

|jy j|| = ||y k || 2 =1), qu’un multiplieur complexe 99 multiplie par la sortie du 

io multiplieur 97 pour adresser la contribution du trajet k a I’additionneur 91 . 

En sortie de I’additionneur 91, les unites 92,93 extraient 

respectivement les parties reelle et imaginaire de y, dont les signes b 1s b 2 
estiment les bits constituant le symbole quaternaire, et dont les valeurs 
absolues foumissent eventuellement des mesures de vraisemblance A 1( A 2 . 

15 Dans le cas d'une modulation BPSK, le module d’estimation de 

symboles 9 peut avoir la meme structure que sur la figure 3, sauf que I’unite 93 
n’est pas necessaire. 

Dans les dessins annexes, le demodulateur est represente sous forme 
synoptique afin de mettre en evidence les differentes etapes de calcul. En 
20 pratique, ce calcul peut etre effectue par un processeur de signal numerique 
programme en consequence, ainsi qu’il est usuel dans les demodulateurs 
numeriques. 

Les inventeurs ont teste les performances du demodulateur ci-dessus 
dans le cas du systeme CDMA retenu pour les reseaux cellulaires de troisieme 
25 generation (voir « UMTS : la generation des mobiles multimedia », P. Blanc et 
a!., L’Echo des Recherches, No. 170, 4 e trimestre 1997 & 1 er trimestre 1998, 
pages 53-68), dans le cas de codes d’etalement de cadence 3,84 Mchip/s avec 
une cadence de symboles de 15 kbit/s (facteur d'etalement SF = 256). Dans le 
systeme considere, le canal de donnees est constitue par la partie reelle (voie 
30 I) d’un signal complexe module en QPSK, transmis par un terminal mobile vers 
une station de base. Le canal pilote correspondant fait partie d’un canal de 
controle constitue par la partie imaginaire (voie Q) du signal complexe module 
en QPSK. Le signal est traite par blocs de 666 ps, a raison de 2560 chips par 
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bloc, ce qui represente 10 bits sur chaque voie avec SF = 256, a savoir 10 bits 
d’information sur la voie I, et 6 bits pilotes et 4 bits de commande sur la voie Q. 
Bien que la modulation radio soit QPSK, les symboles binaires transmis sur 
chaque voie I, Q doivent etre consideres comme faisant I’objet d’une 
5 modulation effective de type BPSK. Dans ces simulations, la longueur de la 
reponse du canal etait d'environ 15 ps (W = 56). Pour chaque bit estime sur la 
voie I a partir d’un signal issu d’une seule antenne, le vecteur r utilise dans la 
relation (1) etait celui resultant de la moyenne des 6 bits pilotes du bloc 
courant, soit (3 = 6. Le canal etait estime au moyen de d = 4 vecteurs propres 
io determines au moyen de M = 3 sous-matrices de dimension S = 6. Pour 
differentes configurations du canal, les simulations ont montre un gain de 1/4 a 
3/4 de decibel sur le rapport Eb/N 0 (Eb etant I'energie par bit) pour un taux 

d’erreur binaire (BER) egal a 10' 2 avant decodage correcteur d’erreurs. Au prix 
d’une augmentation moderee de la complexity du recepteur, ce gain est tres 
appreciable pour augmenter la portee radio du systeme. 



15 
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REVENDICATIONS 

^ Dispositif d estimation de parametres representant un canal de 

propagation radio, comprenant un filtre adapte (5) pour produire des 
observations d un vecteur de reference (r) de W echantillons successifs d’un 
signal forme par filtrage adapte d’un signal (Y) regu selon le canal de 
propagation avec un code d’etalement pilote (Cr), et des moyens (8) 
d’estimation de parametres representant le canal de propagation, incluant des 
vecteurs propres (v k ) d'une matrice d’autocorrelation du vecteur de reference. 

2. Dispositif selon la revendication 1, dans lequel les parametres 
representant le canal de propagation, estimes par les moyens d'estimation de 
parametres (8), incluent en outre des valeurs propres (X k ) de la matrice 

d’autocorrelation du vecteur de reference (r), respectivement associees auxdits 
vecteurs propres (v k ). 

3. Dispositif selon la revendication 1 ou 2, dans lequel les moyens 
d’estimation de parametres (8) comprennent : 

des moyens (80-83) d extraction de M sous-vecteurs de reference (r') 
composes chacun de S echantillons successifs dudit signal forme par 
filtrage adapte, M et S etant des nombres plus grands que 1, les M sous- 
vecteurs de reference correspondant en moyenne a M plages 
energetiques du vecteur de reference (r) ; 

- des moyens (83) d’estimation des matrices d’autocorrelation respectives 

( R ) des M sous-vecteurs de reference ; 

- des moyens (84-85) de determination de sous-vecteurs propres (Vj') 

associes chacun a une valeur propre (\j) et a I’un des M sous-vecteurs 

de reference (r 1 ), chaque sous-vecteur propre etant, pour la valeur propre 
associee, un vecteur propre de la matrice d’autocorrelation estimee du 
sous-vecteur de reference associe ; 

des moyens (85) de formation de chaque vecteur propre estime (v k ) de la 
matrice d’autocorrelation du vecteur de reference en plagant un des 
sous-vecteurs propres (vj 1 ) aux positions des S echantillons du sous- 
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vecteur de reference associe et des echantillons nuls aux W-S autres 
positions. 

4. Dispositif selon la revendication 3, dans lequel les moyens de 
determination des sous-vecteurs propres (85) choisissent lesdits sous-vecteurs 
propres (vj‘), parmi I’ensemble des vecteurs propres des matrices 

f « j 

d’autocorrelation (R ) des M sous-vecteurs de reference (r'), comme etant ceux 

qui sont respectivement associes aux valeurs propres (Xj') de plus grands 
modules. 

5. Dispositif selon la revendication 4, dans lequel les moyens de 
determination des sous-vecteurs propres (85) determinent un nombre donne 
(d) de sous-vecteurs propres (Vj'), associes aux d valeurs propres (Xj‘) de plus 
grands modules. 

6. Dispositif selon I’une quelconque des revendications 1 a 5, dans 
lequel ledit signal regu selon le canal de propagation est compose de plusieurs 
signaux captes au moyen d’antennes differentes. 

7 Demodulateur pour une reception selon un canal de propagation 

radio sur lequel sont definis plusieurs canaux logiques multiplexes incluant un 
canal pilote et au moins un canal de donnees, comprenant des moyens (5) 
d’obtention de plusieurs observations d’un premier vecteur (r) de W 
echantillons successifs d’un premier signal forme par filtrage adapte d'un signal 
(Y) regu selon le canal de propagation au moyen d’un premier code 
d’etalement (Cr) relatif au canal pilote, des moyens (6) d’obtention d’un second 
vecteur (x) de W echantillons successifs d’un second signal forme par filtrage 
adapte dudit signal regu au moyen d’un second code d’etalement (Cx) relatif au 
canal de donnees, des moyens (8) d’estimation de parametres representant le 
canal de propagation, incluant des vecteurs propres (v k ) d'une matrice 

d’autocorrelation du premier vecteur, et des moyens (9) d’estimation de 
symboles pour estimer au moins un symbole d’information porte par le second 
signal sur la base des premier et second vecteurs et desdits parametres 
representant le canal de propagation. 
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8. Demodulateur selon la revendication 7, dans lequel les moyens 
d’estimation de parametres (8) comprennent : 

- des moyens (80-83) d’extraction de M sous-vecteurs de reference (r‘) 
composes chacun de S echantillons successifs du premier signal, M et S 
etant des nombres plus grands que 1 , les M sous-vecteurs de reference 
correspondant en moyenne a M plages energetiques du premier vecteur 
(£); 

- des moyens (83) d’estimation des matrices d’autocorrelation respectives 
(R 1 ) des M sous-vecteurs de reference ; 

- des moyens (84-85) de determination de sous-vecteurs propres (vj') 

associes chacun a une valeur propre (Xj 1 ) et a I’un des M sous-vecteurs 

de reference (r 1 ), chaque sous-vecteur propre etant, pour la valeur propre 
associee, un vecteur propre de la matrice d’autocorrelation estimee du 
sous-vecteur de reference associe ; 

- des moyens de formation de chaque vecteur propre estime (v k ) de la 
matrice d'autocorrelation du premier vecteur en plagant un des sous- 
vecteurs propres (Vj 1 ) aux positions des S echantillons du sous-vecteur de 
reference associe et des echantillons nuls aux W-S autres positions. 

9. Demodulateur selon la revendication 8, dans lequel les moyens de 
determination des sous-vecteurs propres (85) choisissent lesdits sous-vecteurs 
propres (Vj 1 ), parmi I’ensemble des vecteurs propres des matrices 

d'autocorrelation (R 1 ) des M sous-vecteurs de reference (r'), comme etant ceux 

qui sont respectivement associes aux valeurs propres (Xj') de plus grands 
modules. 

10. Demodulateur selon la revendication 9, dans lequel les moyens (85) 
de determination des sous-vecteurs propres determinent un nombre donne (d) 
de sous-vecteurs propres (Vj 1 ), associes aux d valeurs propres (Xj ') de plus 
grands modules. 

11. Demodulateur selon I’une quelconque des revendications 7 a 10, 
dans lequel les parametres representant le canal de propagation, estimes par 
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les moyens d’estimation de parametres (8), incluent en outre des valeurs 
propres (\ k ) de la matrice d’autocorrelation normalisee du premier vecteur, 
respectivement associees auxdits vecteurs propres (v k ). 

12. Demodulateur selon la revendication 11, dans lequel le second 

5 signal porte des symboles d’information binaires (b) modules selon une 
modulation binaire a deplacement de phase, et dans lequel les moyens 
d’estimation de symboles (9) sont agences pour estimer au moins un symbole 
d’information binaire porte par le second signal par le signe (b) de la partie 




10 de vecteurs propres estimes de la matrice d’autocorrelation du premier vecteur, 
v k est un vecteur propre de la matrice d'autocorrelation normalisee du premier 

vecteur, associe a une valeur propre X k , N 0 est une estimation de la puissance 
du bruit sur le canal de propagation, r et x sont respectivement les premier et 
second vecteurs, p est un coefficient predetermine, () H designe le transpose 
15 conjugue, et J|.|| la norme d’un vecteur. 



13. Demodulateur selon la revendication 12, comprenant des moyens 

de determination d’une vraisemblance (A) dudit symbole d’information binaire 
estime (b), proportionnellement a la valeur absolue de la partie reelle dudit 
nombre complexe. 



20 



25 



14. Demodulateur selon la revendication 11, dans lequel le second 

signal porte des symboles d’information quaternaires (b) modules selon une 
modulation quaternaire a deplacement de phase, et dans lequel les moyens 
d’estimation de symboles (9) sont agences pour estimer au moins un symbole 
d’information quaternaire porte par le second signal par deux bits (b 1t b 2 ) 
respectivement definis par les signes des parties reelle et imaginaire du 




ou d est le nombre de 



vecteurs propres estimes de la matrice d’autocorrelation du premier vecteur, v k 



est un vecteur propre de la matrice d'autocorrelation normalisee du premier 
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vecteur, associe a une valeur propre k k , N 0 est une estimation de la puissance 
du bruit sur le canal de propagation, r et x sont respectivement les premier et 
second vecteurs, p est un coefficient predetermine, (,) H designe le transpose 
conjugue, et ||.|| la norme d’un vecteur. 

5 15. Demodulateur selon la revendication 14, comprenant des moyens 

de determination de vraisemblances respectives ( Ai , A 2 ) desdits bits 

( t»i , b 2 ), proportionnellement aux valeurs absolues des parties reelle et 
imaginaire dudit nombre complexe. 

16. Demodulateur selon I’une quelconque des revendications 7 a 15, 

10 dans lequel ledit signal regu selon le canal de propagation est compose de 
plusieurs signaux captes au moyen d’antennes differentes. 





FIG. 4. 
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